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Using experimental data from three different rogue wave supporting systems, determinism and pre¬ 
dictability of the underlying dynamics are evaluated with methods of nonlinear time series analysis. 
We included original records from the Draupner platform in the North Sea as well as time series 
from two optical systems in our analysis. One of the latter was measured in the infrared tail of 
optical fiber supercontinua, the other in the fluence profiles of multifilaments. All three data sets 
exhibit extreme-value statistics and exceed the significant wave height in the respective system by 
a factor larger than two. Nonlinear time series analysis indicates a different degree of determinism 
in the systems. The optical fiber scenario is found to be driven by quantum noise whereas rogue 
waves emerge as a consequence of turbulence in the others. With the large number of rogue events 
observed in the multifilament system, we can systematically explore the predictability of such events 
in a turbulent system. We observe that rogue events do not necessarily appear without a warning, 
but are often preceded by a short phase of relative order. This surprising finding sheds some new 
light on the fascinating phenomenon of rogue waves. 

PACS numbers: 05.45.-a, 47.27.Sd, 42.65.Sf 


Ocean rogue waves are rare events with extreme mag¬ 
nitude and can pose a serious threat to large ships in 
the ocean. Their existence has been questioned for cen¬ 
turies until one such event was actually observed on the 
Draupner platform in the North Sea on January 1, 1995 
DEI- About ten years later, a similar statistical anomaly 
was observed in optical fiber supercontinua (3j. Subse¬ 
quently, the topic of rogue waves received rapidly grow¬ 
ing attention, and heavy-tailed probability density func¬ 
tions characteristic for rogue waves have been reported 
in a multitude of completely different physical systems. 
There appears to be consent that one qualifying crite¬ 
rion for rogue wave statistics is the appearance of large- 
amplitude events at far higher frequency than expected 
by Gaussian statistics [1]. In turn, probability density 
functions exhibit a heavy tail, with a non-vanishing prob¬ 
ability of events with extreme magnitude. A second fre¬ 
quently stated property of rogue waves is their surprising 
appearance and disappearance. Rogue waves have been 
claimed to “appear from nowhere and disappear without 
a trace” [5]. This pictorial statement touches the ques¬ 
tion whether, at least in principle, a prediction of rogue 
waves is possible [3 0 . This discussion has been strongly 
inspired by certain non-stationary solutions of the un¬ 
derlying wave equations, namely Akhmediev breathers 
and the Peregrine soliton isffla, which deterministically 
evolve from a low-level waveform to a transient spike with 
high peak amplitude. Here we tackle the problem of pre¬ 


dictability from a different perspective. We employ non¬ 
linear time-series analysis mm to directly gauge the 
predictability in experimental data. To this end, we use 
a proven method for the detection of even faint traces of 
determinism m- This method enables a quantification 
of rogue wave predictability for the first time. 

We base our analysis on experimental data from three 
different sources. All data sets exhibit a heavytail sta¬ 
tistical distribution of events. The first time series has 
been recorded in the spatial profile of a laser beam that 
underwent multiple filamentation in a gas cell [13] . The 
underlying dynamics is caused by turbulence inside the 
gas, which is induced by local heating due to the plasma 
strings inside the multifilament [Til E> • Time series of 
the spatial fluence distribution were measured with a 
camera in the far field of the multiple filaments. His¬ 
togram analysis of individual camera pixels indicates a 
sub-exponential heavy tail, and common qualifiers for 
rogue wave behavior are exceeded by a large number. 
In this system, we can find events that exceed the sig¬ 
nificant wave height (SWH) by a factor up to ten. The 
SWH is defined as the average of the largest third of 
waves in a record, and surface water waves are usually 
considered rogue when they exceed the SWH by a fac¬ 
tor of ~ 2. The filament data was discretely sampled 
at an acquisition rate of 1 kHz synchronized to the laser 
repetition rate. The data exhibits a linear correlation 
time constant of about 10 ms, which corresponds to ten 
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laser shots. This non-zero correlation time sets the fil¬ 
ament at ion data apart from the behavior of the second 
system, i.e., real-time spectra that have been recorded 
at the output of a fiber supercontinuum m- The latter 
data has been recorded at the much higher megahertz 
repetition rate of a mode-locked oscillator. These spec¬ 
tra exhibit a number of rogue events in their infrared 
tail. We find events that exceed the SWH by a factor of 
4, which clearly qualifies the observed dynamics in the 
fiber as rogue. Finally, the third data set constitutes of 
original surface elevation time series recorded in January 
1995 on the Draupner oil platform mm- These rela¬ 
tively short data sets contains rogue events that exceed 
the SWH by a factor of 2.15 to 2.3. 

Table [I] shows a comparison of the different parameters 
of the three data sets. We include the shape parameter 
/3 of the histograms of the respective physical quantity x, 
as determined by a fit of the probability density function 
(PDF) to the Weibull function jTf] 
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Both optical system exhibit sub-exponential PDFs with 
minimum values of /3 < 1, while the ocean records exhibit 
/3 ss 2. All distributions are far away from normally dis¬ 
tributed data, which would be indicated by /3 > 10. The 
dynamics in all systems exceed the SWH or an equiva¬ 
lently defined measure by a factor >2. In this linear sta¬ 
tistical comparison of the data, there is another slightly 
more hidden aspect that sets the fiber data apart from 
the other systems. The temporal autocorrelation of the 
fiber system is essentially a delta function, whereas the 
other two time series exhibit autocorrelation widths i corr 
that are clearly larger than the sampling interval St. Fig¬ 
ure [T] shows characteristic rogue events for the three sys¬ 
tems as well as their autocorrelation functions. To sim¬ 
plify inter-system comparison, time axes in Figs. 3a, b) 
have been scaled in units of the autocorrelation width. 


Data set 

N 

St 

Tcorr 

P 

Excess SWH 

Fiber 

5,000 

40 ns 

<1 ps 

0.87 

4.2 

Filament 

60,000 

1 ms 

12 ms 

0.44 

10 

Draupnerl 

2,520 

0.47s 

20 s 

2.04 

2.15 

Draupner2 

2,520 

0.47s 

20 s 

1.8 

2.3 


TABLE I: Intersystem comparison of rogue wave characteris¬ 
tics. N : total number of samples in time series. 5t : sampling 
interval. t c orr : linear autocorrelation time (1/e width). /3 : 
Shape factor determined by Weibull fit to the histograms [ 13 - 
Excess SWH: Maximum wave height observed relative to sig¬ 
nificant wave height. Draupnerl was recorded on 01/01/1995 
starting at 3:20pm; Draupner2 on 01/19/1995 at 11pm. 

For a deeper analysis of the time series records, we 
employed one of the standard methods of nonlinear 
time series analysis, namely the Grassberger-Procaccia 
algorithm (GPA) [IT]. Let us suppose that a physi¬ 
cal quantity x has been sampled at a sampling inter- 



FIG. 1: Examples of rogue events (a) in multifilament dy¬ 
namics Q3], (b) the oceanographic context (Draupnerl data 
set [L m ) i and ( c ) soliton dynamics in a nonlinear fiber [3]. 
Abscissae in (a) and (b) have been rescaled to equal units 
in terms of autocorrelation width, see insets. Solid lines: 
long-term averages. Dashed horizontal lines: significant wave 
height (SWH). Whenever applicable, the linear 1/e autocor¬ 
relation width is shown as dashed vertical lines in the insets. 

val St with a record length N. This yields a time se¬ 
ries x that we denote by {aq, aq, x 3 , ...,x^} in the fol¬ 
lowing. From these time series, we select sub-series 
y i, m = {xi,x i+ i,...,x i+m }. Here the record length m 
is called the embedding dimension. These sub-series are 
then analyzed according to their Euclidian distances 
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which are accumulated for all i and j > i in a histogram 


C m (r) = 
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where e determines the bin size. Deterministic behavior is 
indicated by anomalously large counts C m (r) at r = 0.01 
to 0 . 1 r max , where r max denotes the largest non-zero bin of 
the histogram. In an illustrative way of speaking, such 
an anomaly can be understood as “deja-vus”, i.e., the 
frequent appearance of certain time sequences y that are 
beyond a completely random distribution of events. De¬ 
tection of such a statistical anomaly may be fairly ob¬ 
vious in physical systems with a low internal number of 
degrees of freedom, see for example [18]. In turbulent 
systems, however, it requires very careful checks to judge 
the presence of determinism in a data set. As a null 
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hypothesis, one therefore compiles surrogate data sets 
that resemble the original data set in all aspects of linear 
descriptive statistics. In particular, these surrogates fea¬ 
ture an identical histogram C m (r) and also exhibit nearly 
the identical power spectrum and correlation function. 
We generated appropriate surrogate data in a three-step 
procedure, starting with the amplitude adjusted Fourier 
transform originally proposed by Theiler et al. 119| . To 
remove any residual bias at low Fourier frequencies, we 
applied an iterative method suggested by Schreiber and 
Schmitz |20| . Finally, we apply an annealing method 
to further improve agreement of the linear correlation 
function between original and surrogate ED- All three 
steps reorder the time series without affecting the his¬ 
togram. The randomization procedure heavily relies on 
the quality of true random number sequences that need 
to be free from any residual correlation. This constraint 
immediately rules out algorithmic random number gen¬ 
erators. We therefore resorted to two approved sources 
of true random data [22J [22], which have been generated 
by hardware generators. These random number genera¬ 
tors are based on Johnson noise of a resistor or similar 
sources of uncorrelated noise. Additionally, both sources 
have been carefully checked to comply with a suite of 
statistical torture tests [ 22 j . 

Figure [2] shows results of our surrogate analysis for 
the three different rogue wave systems. As optical rogue 
waves only appear in a fraction of the spectrum or beam 
profile, we have selected the respective time series that 
exhibit the most pronounced heavy tail, as indicated by 
the smallest value of the shape parameter /3. Figure [2] 
shows exemplary histograms C m (r) for the original data 
sets as well as for 100 surrogates. In these examples, the 
record length was adapted to the equivalent of one corre¬ 
lation length; for the fiber soliton system with its vanish¬ 
ing correlation length, we employed m = 4. It needs to 
be emphasized that the below conclusions are reproduced 
for a large range of record lengths. Independent of m and 
the system under consideration, there are no statistically 
significant deviations between original and surrogates at 
large r as expected. Determinism in the data shows up 
as a lower histogram count C m (r) for the surrogates at 
low r. This behavior appears most pronounced for the 
multifilament system [Fig. [2^a)] , with a statistical signif¬ 
icance reaching values of 100 cr surr well beyond the auto¬ 
correlation time of the system (equivalent to m = 12 ). 
This finding further supports small-scale turbulence as 
the origin of the observed dynamics [MIS]. In contrast, 
the fiber optical system [Fig. [2jc)] exhibits maximum sta¬ 
tistical significances in the vicinity of 3cr SU r r . At larger 
values of m and at low values of r, the surrogate his¬ 
tograms tend to lie even above the original histogram, 
which has to be interpreted as stochastic behavior, i.e., 
an absence of determinism in the data. It therefore seems 
to be evident that the fiber system is driven by a com¬ 
pletely uncorrelated noise source. Considering the under¬ 
lying physics of the system, amplified spontaneous emis¬ 
sion noise of the pump laser [25, [ 2 B] appears the only 
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FIG. 2: Nonlinear correlation analysis of rogue-wave systems. 
Shown is the histogram count C(r) as a function of Euclidian 
distance r, cf. Eq. Red line: histogram of original data 
set. Blue line: average of surrogates. Dashed blue lines: sur¬ 
rogate standard deviation spread. Embedding dimensions m 
have been adapted to the correlation widths in (a) the nrulti- 
filament system and (b) an analysis of the Draupnerl event, 
(c) Analysis for solitons in a nonlinear fiber, m = 4 has been 
chosen to avoid trivial behavior, (d) Maximum detected sta¬ 
tistical significance of deterministic behavior as a function of 
m. Symbols indicate the position of the histograms (a-c) in 
the plot. Significances < 3 can safely be ignored as they 
indicate stochastic behavior within the expected credibility 
limits. Green curve: multifilament. Dark blue curve: Draup¬ 
nerl data set. Light blue curve: Draupner2 data set. Red 
curve: fiber solitons. 

possible candidate for such a noise source. The nonlinear 
fiber propagation itself is a completely deterministic pro¬ 
cess, which is ruled by the nonlinear Schrodinger equa¬ 
tion. Nevertheless, this system is infamous for its large 
noise amplification [251 EZ|, which gives rise to a rapid 
loss of coherence during the highly nonlinear propaga¬ 
tion. Analysis of the Draupner data sets [Fig. [2)(b)] is in¬ 
dicative of some determinism in the data. However, this 
determinism does not extend much beyond T CO rr (equiv¬ 
alent to m ~ 40). Such a weak deterministic character 
of ocean surface waves has been reported before [24] , yet 
without including rogue waves in the analysis. This be¬ 
havior is certainly different from the stochastic behavior 
in Fig. [2j)c) . As an important conclusion, we find that 
the degree of determinism is completely independent of 
the extreme-value statistics, i.e., there is no correlation 
between the parameter [3 and the predictability in the 
system. Moreover, there actually seem to be at least two 
different classes of rogue waves, those that emerge dur¬ 
ing turbulence and a second class that is generated by 
quantum noise. 
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FIG. 3: Histogram GPA analysis of time series segments in 
immediate vicinity of a rogue event in the multifilament data 
(red curve) in comparison to arbitrarily chosen segments of 
equal length (solid blue curve), (a) Pre-rogue statistics. Ref¬ 
erence time series immediately preceding 289 isolated rogue 
events, (b) Post-rogue statistics. The statistical spread for 
arbitrarily chosen data (100 data sets) is indicated by the 
dashed blue curves. 

So far, we have drawn conclusions from time series that 
only contained a relatively small fraction of rogue events. 
In principle, it is conceivable that these rogue events are 
generated by a second underlying mechanism in the sys¬ 
tem, which does not exhibit the same determinism as 
the rest of the trace. Given the very few records of real 
ocean rogue waves, this issue appears difficult to decide, 
despite some attempts towards a systematic collection of 
ocean rogue events [28] . In contrast, the large number 
of rogue events in our turbulent multifilament system 
provides a unique opportunity to tackle this question. 
To isolate the dynamics accompanying rogue events, we 
select sub-series m and y2 j°i ^ immediately pre¬ 

ceding or following a rogue event, respectively. The index 
i was chosen such that it either marks the beginning or 
the end of a single rogue event exceeding the SWH by a 
factor 2. We end up with 289 carefully chosen rogue 
events in our data, ensuring that they occurred more 
than one second and 260 fim on the camera chip apart 
from each other. For the null hypothesis, we selected 
time series for randomly chosen i at the same pixels of 
the camera. We then independently compiled histograms 
of the rogue-wave related sub-series and compared them 
to the histograms of the random data sets. The results 
are shown in Fig. [3] Looking again at the low-r part of 
the histograms at about 20% of r max , we notice a statis¬ 
tically significant anomaly similar to the one in Fig. [2] 
i.e., at certain values of r we notice up to 5 times as 
many counts C m (r) for the rogue data sets than for arbi¬ 
trarily chosen data. Nearly the same statistical anomaly 
is observed for the segments right after a rogue event 
[Fig. [3](b)]. This means that the appearance of a rogue 
event is accompanied by characteristic dynamics in the 
time series. Similar peculiarities were observed in a wide 
range of fields, reaching from the appearance of epileptic 
seizures in medicine [ 55 ] to geomagnetic storms [ 30 ] . The 


concomitant contraction of the parameter space can be 
taken as an analogy to the famous calm that precedes a 
storm. In nonlinear dynamics, such phenomena are typi¬ 
cally understood as a shrinking of phase space connected 
to the appearance of the event, which, in turn, indicates 
the existence of at least one strange attractor. Most im¬ 
portantly, this finding appears to contrast the statement 
that “rogue waves appear from nowhere and disappear 
without a trace”. 

The observed predictability is rather robust, even if we 
select shorter sub-series at one or two correlation times 
distance to the rogue event. While the statistical signif¬ 
icance certainly drops compared to Fig. [3] we only start 
to lose evidence for predictability when we select seg¬ 
ments more than ss 5r CO rr away from the rogue event. 
One can certainly not directly transfer this finding to 
ocean rogue waves. At best one may expect to predict an 
ocean rogue wave a few ten seconds before impact, and it 
would require many future sightings to isolate character¬ 
istic patterns preceding an ocean rogue wave. Therefore 
any practical rogue wave prediction appears not overly 
realistic, despite the determinism in the system. 

At first glance, the complete absence of determinism in 
one of the systems may actually be surprising as all three 
considered systems are described by variants of the Non¬ 
linear Schrodinger Equation (NLSE) and exhibit similar 
soliton-like solutions. This formal similarity has always 
been considered a unifying aspect of rogue-wave support¬ 
ing systems ED- The appearance of turbulence presup¬ 
poses a spatially extended character of the system, and 
it is striking that both the ocean and the multifilament 
system are two-dimensional systems whereas the nonlin¬ 
ear fiber scenario is a one-dimensional system, which can 
also show turbulent behavior [321 EE]- Nevertheless, the 
driver for rogue waves in the fiber system is quantum 
noise of the supercontinuum pump laser, giving rise to a 
complete lack of determinism in the observed dynamics. 
In some sense of speaking, therefore, extreme events in 
the fiber optic context appear as the better rogue waves 
because of their guaranteed unpredictability. In the orig¬ 
inal ocean context, in principle, rogue waves seem to be 
predictable on short time scales, which further demysti¬ 
fies them as a deterministic phenomenon. 

We gratefully acknowledge Daniel Solli (UCLA, Los 
Angeles, California) for providing original rogue wave 
data for the fiber-optical system and Statoil ASA for the 
permission to use the data recorded on the Draupner oil 
platform in January 1995. Mads Haahr is acknowledged 
for giving us generous access to true random number se¬ 
ries at random.org. Moreover, we thank Michel Olagnon 
(IFREMER, Brest, France), Sverre Haver (Norwegian 
University of Science and Technology, Trondheim), and 
Anne Karin Magnusson (Norwegian Meteorological In¬ 
stitute, Oslo) for their support. We further acknowledge 
Fedor Mitschke (University of Rostock) and Claus Rop¬ 
ers (University of Gottingen) for helpful discussions. 









5 


[1] S. Haver and D. Karunakaran, “Probabilistic description 
of crest heights of ocean waves,” in: Proc. 5 th Interna¬ 
tional Workshop of Wave Hindcasting and Forecasting, 
Melbourne, FL (1998). 

[2] S. Haver, “A possible freak wave event measured at the 
Draupner Jacket January 1 1995,” in: Proc. of Rogue 
Waves, M. Olagnon and M. Prevosto (eds.), Brest, France 
(2004). 

[3] D. R. Solli, C. Ropers, P. Koonath, and B. Jalali, “Op¬ 
tical rogue waves,” Nature 450 , 1054-1057 (2007). 

[4] N. Akhmediev and E. Pelinovsky, “Editorial Rogue 
Waves Towards a Unifying Concept?: Discussions and 
Debates,” Eur. Pliys. J. Special Topics 185, 1-4 (2010). 

[5] N. Akhmediev, A. Ankiewicz, and M. Taki, “Waves that 
appear from nowhere and disappear without a trace,” 
Phys. Lett. A 373 , 675-678 (2009). 

[6] M.-R. Alam, “Predictability horizon of oceanic rogue 
waves,” Geophys. Res. Lett. 41 , 8477-8485, (2014). 

[7] M. Olagnon and M. Prevosto, “Are rogue waves be¬ 
yond conventional predictions?” Proc. 14 th ’Aha Hu- 
liko’a Hawaiian Winter Workshop (2005). 

[8] B. Kibler, J. Fatome, C. Finot, G. Millot, F. Dias, 

G. Genty, N. Akhmediev, and J. M. Dudley, “The Pere¬ 
grine soliton in nonlinear fibre optics,” Nature Phys. 6, 
790-795 (2010). 

[9] N. Akhmediev, A. Ankiewicz, J.M. Soto-Crespo, 
J.M. Dudley, “Rogue wave early warning through spec¬ 
tral measurements?” Phys. Lett. A 375 , 541-544 (2011). 

[10] A. Chabchoub, N. Hoffmann, H. Branger, C. Kharif, 
and N. Akhmediev, “Experiments on wind-perturbed 
rogue wave hydrodynamics using the Peregrine breather 
model,” Phys. Fluids 25 , 101704 (2013). 

[11] P. Grassberger and I. Procaccia, “Characterization of 
strange attractors,” Phys. Rev. Lett. 50 , 346-349 (1983). 

[12] H. D. I. Abarbanel, R. Brown, J. J. Sidorowich, and 
L. S. Tsimring, “The analysis of observed chaotic data 
in physical systems,” Rev. Mod. Phys. 65, 1331-1392 
(1993). 

[13] S. Birkholz, E. T. ,J. Nibbering, C. Bree, S. Skupin, 
A. Demircan, G. Genty, and G. Steinmeyer, “Spatiotem- 
poral Rogue Events in Optical Multiple Filamentation,” 
Phys. Rev. Lett. Ill, 243903 (2013). 

[14] Y.-H. Cheng, J. K. Wahlstrand, N. Jhajj, and 

H. M. Milchberg, “The effect of long timescale gas dy¬ 
namics on femtosecond filamentation,” Opt. Express 21, 
4740-4751 (2013). 

[15] A. Ryabtsev, S. Pouya, M. Koochesfahani, and M. Dan- 
tus, “Vortices in the wake of a femtosecond laser fila¬ 
ment,” Opt. Express 22 , 26098-26102 (2014). 

[16] D.R. Solli, G. Herink, B. Jalali, and C. Ropers, “Fluctua¬ 
tions and correlations in modulation instability,” Nature 
Photon. 6, 463-468 (2012). 

[17] W. Weibull, “A statistical distribution function of wide 
applicability,” J. Appl. Mech. Trans. ASME 18, 293-297 
(1951). 


[18] E. N. Lorenz, “Deterministic Nonperiodic Flow,” J. At¬ 
mos. Sci. 20 , 130-141 (1963). 

[19] J. Theiler, S. Eubank, A. Longtin, B. Galdrikian, and 
J. D. Farmer, “Testing for nonlinearity in time series: the 
method of surrogate data,” Physica D 58, 77-94 (1992). 

[20] T. Schreiber and A. Schmitz, “Improved Surrogate Data 
for Nonlinearity Tests,” Phys. Rev. Lett. 77 , 635-638 
(1996). 

[21] T. Schreiber, “Constrained randomization of time series 
data,” Phys. Rev. Lett. 80 , 2105-2108 (1998). 

[22] G. Marsaglia, “The Marsaglia Random Num¬ 
ber CD ROM including the Diehard Battery of 
Tests of Randomness,” (1995), downloaded from 
http: //stat.fsu.edu/pub/diehard/ 

[23] M. Haahr, True Random Number Service, 
https://www.random.org, surrogates have been based 
on the January 2014 data sets on this server. 

[24] S. Elgar and G. Mayer-Kress, “Observations of the fractal 
dimension of deep- and shallow-water waves,” Physica D 
37 , 104-108 (1989). 

[25] K. L. Corwin, N. R. Newbury, J. M. Dudley, S. Coen, 
S. A. Diddams, K. Weber, and R. S. Windeler, “Funda¬ 
mental Noise Limitations to Supercontinuum Generation 
in Microstructure Fiber,” Phys. Rev. Lett. 90 . 113904 
(2003). 

[26] J. M. Dudley, G. Genty, and S. Coen, “Su¬ 
percontinuum generation in photonic crystal fiber,” 
Rev. Mod. Phys. 78 , 1135-1183 (2006). 

[27] N.R. Newbury, B.R. Washburn, K.L. Corwin, and 

R. S. Windeler “Noise amplification during supercontin¬ 
uum generation in microstructure fiber,” Opt. Lett. 28 , 
944-946 (2003). 

[28] M. Christou and K. Ewans, “Examining a comprehen¬ 
sive data set containing thousands of freak wave events.” 
Proc. ASME 2011 30th Int. Conf. Ocean Offshore Arctic 
Engin. (OMAE2011). 

[29] J. Martinerie, C. Adam, M. Le Van Quyen, M. Baulac, 

S. Clemenceau, B. Renault, and F.J. Varela., “Epileptic 
seizures can be anticipated by non-linear analysis” Na¬ 
ture Medicine 4 , 1173 (1998). 

[30] C.-C. Chen, Y.-J. Chuo, F.-L. Wang, H.-Y. Yen, and 
C.-H. Chen, “Correlation Dimension and Its Temporal 
Variations in Geomagnetic Total Field During Storms,” 
Terr. Atmos. Ocean Sci. 16 , 435-443 (2005). 

[31] J. M. Dudley, F. Dias, M. Erkintalo, and G. Genty, “In¬ 
stabilities, breathers and rogue waves in optics,” Nature 
Photon. 8, 755-764 (2014). 

[32] F. Mitschke, G. Steinmeyer, and A. Schwache, ’’Gener¬ 
ation of one-dimensional optical turbulence,” Physica D 
96 , 251-258 (1996). 

[33] V. Zakharov, F. Dias, and A. Pushkarev, “One¬ 
dimensional wave turbulence,” Phys. Rep. 398, 1-65 
(2004). 




